arXiv:1506.03646vl [physics.atom-ph] llJun2015 


Influence of the Coulomb potential on above-threshold ionization: a quantum-orbit 
analysis beyond the strong-fleld approximation 

X.-Y. C. Poli^, H. Schomerus^ and C. Figueira de Morisson Faria^ 

^Department of Physics and Astronomy, University College London, 

Gower Street, London WCIE 6BT, United Kingdom 
^ State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 

Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China 
® Department of Physics, Lancaster University, LAI 4YB Lancaster, United Kingdom 

(Dated: June 12, 2015) 

We perform a detailed analysis of how the interplay between the residual binding potential and a 
strong laser field influences above-threshold ionization (ATI), employing a semi-analytical. Coulomb- 
corrected strong-held approximation (SFA) in which the Coulomb potential is incorporated in the 
electron propagation in the continuum. We Hnd that the Coulomb interaction lifts the degeneracy of 
some SFA trajectories, and we identify a set of orbits which, for high enough photoelectron energies, 
may be associated with rescattering. Furthermore, by performing a direct comparison with the 
standard SFA, we show that several features in the ATI spectra can be traced back to the inhuence 
of the Coulomb potential on different electron trajectories. These features include a decrease in the 
contrast, a shift towards lower energies in the interference substructure, and an overall increase in 
the photoelectron yield. All features encountered exhibit a very good agreement with the ab initio 
solution of the time-dependent Schrodinger equation. 
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I. INTRODUCTION 

When matter interacts with a strong laser field of peak 
intensity around 10^"^ W/cm^, the outmost electron may 
be freed by absorbing many more photons than necessary. 
This very highly nonlinear process is known as above¬ 
threshold ionization (ATI) and has attracted consider¬ 
able attention since the early work of Agostini and co¬ 
workers [I|; for a review see Ref. 0. For typical param¬ 
eters employed in experiments, i.e., near infra-red laser 
fields, it is commonly accepted that the electron reaches 
the continuum by tunnel ionization. If the released elec¬ 
tron revisits the parent ion in the presence of the laser 
field di , this results in various additional highly nonlin¬ 
ear phenomena, such as high-order ATI (HATI) Q : high- 
order harmonic generation (HHG) Q, and nonsequen¬ 
tial double ionization (NSDI) Q. Recently, considerable 
progress has been made in the study of these nonlin¬ 
ear strong-field phenomena. For example, both ATI and 
HHG have been employed as an important technique to 
explore the electron shell structure and sub-femtosecond 
dynamics [sl-fl^. and NDSI has opened the door to the 
study of strong-field electron-electron correlation [l3l - [l^ . 

In order to uncover the underlying physics of these 
highly nonlinear phenomena, many theories and mod¬ 
els have been proposed, such as the ab initio solution 
of the time-dependent Schrodinger equation (TDSE) 0 
and the quantum-orbit theory within the strong-field ap¬ 
proximation (SFA) [Tsl - I^ . Since the TDSE contains no 
physical approximation, its outcome is widely taken as a 
benchmark to evaluate the data in experiments and the 
calculations of other theories and models [2ll - [^ . How¬ 
ever, in many occasions the TDSE does not provide a 
transparent physical picture. Furthermore, since the nu¬ 
merical effort involved in ab initio computations increases 


exponentially with the degrees of freedom, its implemen¬ 
tation is impractical for strongly correlated multielectron 
systems. In contrast, the quantum-orbit theory provides 
very clear physical insight in terms of distinct electron 
ionization trajectories, and its outcome is qualitatively 
consistent with the experimental data. Therefore, it has 
been widely and successfully used in the modeling of 
strong-field phenomena |24 - I^ 

One should note, however, that the validity of the 
conventional quantum-orbit theory is limited. In fact, 
the use of the SFA before the application of the saddle- 
point approximation implies that a considerable amount 
of physics is left out for the sake of a clear and intu¬ 
itive picture 0. In particular, the SFA fully neglects the 
effect of the Goulomb potential of the parent ion on the 
ionized electrons, approximating the continuum states by 
field-dressed plane waves For single charged nega¬ 
tive ions, this approximation is justified as the Goulomb 
interaction between the neutral core and the freed elec¬ 
tron is absent. However, for atoms and molecules, this 
interaction is present, so that the SFA only works qual¬ 
itatively. Furthermore, in recent years, several features 
have been observed which clearly highlight the influence 
of the Coulomb potential. Examples are the so-called 
low-energy structure (LES) in ATI spectra [23 - 1^ . fan¬ 
shaped structures in photoelectron momentum distribu¬ 
tions [Sll - fs^ , and the violation of the fourfold symmetry 
in angular electron distributions for elliptically polarized 
fields [3^[^. 

Motivated by these observations, many methods have 
been developed in the past few years in order to account 
for the Coulomb potential in orbit-based methods. These 
include (i) using Coulomb-Volkov functions to describe 
the electron continuum states in the SFA [38l - l40l |: (ii) 
incorporating the binding potential in the electron prop- 
agation using the eikonal Volkov approximation |4ll . l42j| ; 
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(iii) a Coulomb-corrected SFA (CCSFA), which takes the 
trajectories from the SFA theory as a zero-order approx¬ 
imation and accounts for the Coulomb field perturba- 
tively [11,141,lU; (iv) a quantum-trajectory Monte Carlo 
(QTMC), which is based on Classical-Trajectory Monte- 
Carlo (CTMC) simulations, but considers the phase of 
each trajectory [45l |: and (v) initial-valu e rep resentations 


such as the Herman-Kluk propa gato r 
Coupled Coherent States method |5Cll452|. 


and the 


Most of the above-mentioned approaches have been ap¬ 
plied to and tested on direct ATL This phenomenon is 
a particularly good testing ground for Coulomb correc¬ 
tions for two main reasons. First, the momentum range 
involved is relatively low, so that the influence of the 
Coulomb potential is expected to be significant. Sec¬ 
ond, in contrast to high-order ATI, hard collisions are ex¬ 
pected to be absent, so that the Coulomb corrections are 
in principle easier to implement. In particular the influ¬ 
ence of the Coulomb potential on quantum-interference 
patterns has attracted a great deal of attention . Fur¬ 
thermore, it has been shown that the presence of the 
Coulomb potential considerably alters the topology of 
the orbits, giving rise to types of trajectories that are 
absent in the SFA [d^. In particular, in Ref. [s^ it has 
been shown that sub-barrier corrections are necessary in 
order to obtain the correct phases in the ATI electron 
momentum distributions. 


In this work, we develop a quantum-orbit theory with 
Coulomb interactions that, besides the effect on the 
phase, also accounts for the influence of the interactions 
on the semiclassical amplitudes. We find that the ampli¬ 
tude is significantly modified both via the atomic dipole 
moment at ionization and due to the altered stability of 
the trajectories during the electron propagation through 
the continuum. This Coulomb-corrected method is then 
employed to study the influence of the Coulomb potential 
on the direct ATI ionization spectrum of Hydrogen. We 
perform a systematic investigation of how the Coulomb 
coupling changes the topology of the trajectories, which 
are either decelerated or accelerated with regard to their 
SFA counterparts. This leads to a decrease in the phase 
difference between the contributions from different types 
of trajectories, which influences the interference patterns 
in the spectra. We also discuss how momentum non¬ 
conservation lifts the degeneracy of certain SFA trajec¬ 
tories. Furthermore, we verify that the distinction be¬ 
tween direct and rescattered trajectories is blurred by 
the presence of the binding potential, which causes a set 
of trajectories to go around the core. 

Our results show that the spectrum calculated with 
this method is in much better agreement with the ab 
initio TDSE result than the predictions of the standard 
SFA. In particular, the Coulomb-corrected theory recov¬ 
ers the much weaker contrast in the interference substruc¬ 
ture observed in the TDSE, and relates this effect to the 
unequal semiclassical weights of the electron trajectories 
in the presence of the Coulomb interactions. Similarly 
to what has been encountered in [35j |. we also observe 
that the positions of the interference maxima in the spec¬ 


trum from the quantum-orbit theory and TDSE result 
are shifted with respect to the SFA simulations. How¬ 
ever, our model indicates that these shifts mainly stem 
from the modified electron propagation in the continuum. 

This article is organized as follows. In Sec. HIl we de¬ 
scribe the theoretical models employed in this paper, 
namely the standard SFA and the Coulomb-corrected 
SFA, starting from the TDSE. In Sec. IHII we apply the 
theories to direct ATI and discuss the consequences of the 
Coulomb interactions, first in terms of the individual tra¬ 
jectories and then for the resulting ionization spectrum. 
Finally, in Sec. lIVl we summarize the main conclusions to 
be drawn from this work. We use atomic units through¬ 
out. 


II. THEORETICAL MODELS 

The underlying framework for the subsequent discus¬ 
sions is the time-dependent Schrodinger equation 

idt\'il}{t)) = . ( 1 ) 

In the ionization problems considered in this work, the 
Hamiltonian separates into two parts, H[t) = Ha+Hj{t). 
Here 


Ha = ^ + V{r) (2) 

denotes the field-free one-electron atomic Hamiltonian 
and the hats denote operators. In the problem addressed 
here, we consider a Coulomb-type potential 


V{r) 



(3) 


where 0 < C < 1 is an effective coupling, which we vary 
in a continuous fashion in order to assess the influence of 
the Coulomb potential. For Hydrogen, (7=1. Further¬ 
more, Hi{t) describes the interaction with the laser field. 
In the velocity and length gauges, this interaction is given 
by Hi{t) = p- A(t)-f A^/2 and Hi{t) = — f •E(t), respec¬ 
tively, where E(t) = —dA{t)/dt is the external laser field. 
The length gauge provides us with the physical picture 
of ionization as a tunneling process driven by an effective 
time-dependent potential. This gauge will be employed 
throughout. 

The time-evolution operator associated with this 
Hamiltonian is of the general form 


U{t,to) = T exp 



H{t')dt' , 


(4) 


where T denotes time-ordering. This operator takes a 
wave function from a time to to a time t, i.e., \ip{t)) = 
[/(t,to)|'0(to)), and satisfies 


idtU{t,to) = H{t)U{t,to), 

-idt,U(tM = U(tM)H{to). (5) 
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Employing the Dyson equation, the time-evolution oper¬ 
ator may be written as 

t 

U{t,to) = Ua{t,to) - i j , (6) 

^0 

where Ua(t,tQ) is the time-evolution operator associated 
with the field-free Hamiltonian. 

For above-threshold ionization, the initial state is a 
bound state |'0o): while the final state is a continuum 
state \tpp{t)) with drift momentum p. This gives the 
ionization amplitude 0 

t 

M(p) =-i lim f dt', 

t—^OO J 

— oo 

(7) 

which is formally exact. 


A. Strong-field approximation 


Equation 0 cannot be solved in closed form, so that 
approximations are required in order to compute the ATI 
transition amplitude via analytical methods. A popular 
approximation is to replace by the Volkov time 

evolution operator in Eq. 0. This implies 

that the continuum has been approximated by Volkov 
states, i.e., by field-dressed plane waves \ip^\t)), where 


(t)) = (r|p(t))exp 






( 8 ) 


with 


(r|p(t)) 


exp [ip (t) • r] 
(27r)3/2 


(9) 


Here p(t) = p-b A(t) in the length gauge and p(<) = p in 
the velocity gauge, so that 

[soil . This is the key idea behind the str ong -field approx¬ 
imation or Keldysh-Faisal-Reiss theory . For 

detailed discussions see, e.g., and the recent tu¬ 

torials ( 5 ^ . Within the SEA, the amplitude (l7|) is then 
given brn, US min 


M{p)=-1 J dt'{p + A{t')\Hj{t')\^o)e^^^^’^'l ( 10 ) 

— 00 

Here 

1 

S{p,t') = --J [p + A{T)]‘^dT + Ipt' (11) 

is the semiclassical action, where Ip gives the ionization 
potential and A{t) denotes the vector potential of the 
laser field. In Eq. m, we have also employed the nota¬ 
tion IV'o(i')) = lipo)- 


For sufficiently high intensity and low frequency of the 
laser field, the temporal integration in Eg. (HM can be 
evaluated by the saddle-point method [13, l23i which 
seeks solutions such that the action m is stationary. 
The corresponding saddle-point equation reads 


[P + A(t')]^ 
2 


+ Ip — 0. 


( 12 ) 


Physically, Eq. m ensures the conservation of energy at 
the ionization time t', which leads to complex solutions 
tg. In terms of these solutions, the transition amplitude 
OT can then be written as 


M(p) ^ ^C(4) (p + Aits) \Hjits)\ V'o) , 


S 


where the prefactors 


(13) 


Cits) 


2'Ki 


d^sip,ts)/dt^s 


(14) 


are expected to vary much more slowly than the action 
for the saddle-point approximation to hold. Since each 
solution ts represents a distinct trajectory of the elec¬ 
tron in the laser field, the sum in Eq. m denotes the 
interference between different quantum paths, which has 
been extensively studied in the literature (for reviews see, 
e.g., UM)- One should note that in the SEA, the field- 
dressed momentum is conserved. 


B. Coulomb-corrected SFA 

Within the SFA, an electron no longer feels the atomic 
potential after it has been promoted into the continuum 
at the time tg, resulting in the considerable deviations 
between this model and experimental results. In this 
section we describe a Coulomb-corrected SFA (CCSFA) 
which cures this shortcoming of the SFA (for similar ap¬ 
proaches see [ 3 ^ Idll - li^ l. 

First we note that in presence of the Coulomb poten¬ 
tial, momentum is no longer conserved, and the time evo¬ 
lution operator depends on both f and p. As a result, 
the time evolution operator cannot be diagonalized by 
the Volkov states ([5]). 

Instead, adopting Feynman’s path integral formalism 
the Coulomb corrected transition amplitude reads 

ft fPit) f Pj. 

M(p) = -i lim / dt' Vp 

t^°oJ-oo Jp(t') J (27r)^ 

xe*^(P''-’‘''‘)(po + A(0 |Fl7(OI V'o) , (15) 

where the action is given by 

-^(P, r, f', t) = Ipt' + Sip, r, t', t), (16) 

with 

Sip, r, t', t) = - [r(T) . p(t) -t- i?(r(r), p(t), r)]dr, 

( 17 ) 
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and 

i7(r(T), p(r), r) = i [p(r) + A(t)]^- ^ 

4 y^r(T) • r(r) 

(18) 

One should note that the problem is solved in the length 
gauge, and Eq. (HI can be obtained from the standard 
length-gauge Hamiltonian by a partial integration. This 
issue has been discussed in detail in Ref. [6^ . 

Following the same procedure as for the SEA, we can 
now obtain the Coulomb-corrected transition amplitude 
by applying the saddle-point approximation. By con¬ 
struction, the saddle-point equation on t' leads to the 
condition 


[P + A(t')]^ 
2 


+ Ip + E[r(t^)] — 0. 


(19) 


The Coulomb-corrected transition amplitude is then 
given by 


M(p) = —i 




Vr 


X {p{t,)+AiU) \Hiits)\ V'o) . 


( 20 ) 


We then use the semi-classical path integral theory de¬ 
veloped by Van Vleck and Gutzwiller on the paths 
[p(t)], and [r(t)]. The associated saddle-point equations 
take the form of classical equations of motion for the tra¬ 
jectories. 


p(r) = -VrU[r(T)], 

(21) 

f(T) = p(t) + A(t) . 

(22) 


The solutions of these equations in general are again com¬ 
plex due to the additional condition (1191) . 

In terms of such solutions, the Coulomb-corrected 
transition amplitude finally reads 


M{p) = 

S 


2TTiy^S"{ts) 


det 


dPsjt) 

drs{ts) 



X (psits) + A(G) V^o) . (23) 


The sum is over the classical trajectories that begin at 
position r{ts) at time tg, and end at momentum p{t) at 
time < —> oo (sum over multiple solutions with identical 
tg is implied). Each trajectory contributes a term with 
phase given by the Maslov index i/ and the classical action 
S{rg,Pg,tg,t) given by Eq. (|T6)l at the stationary values 
ts,pg,rs. In practice, t is defined at the end of the pulse, 
which should be taken to be sufficiently long. 

Due to the presence of the binding potential, the 
above-stated equation exhibits branch cuts for Re[r 5 (r) • 
Ps(t)] < 0 and Im[rs(r) • rs(T)] = 0. For vanishing trans¬ 
verse momenta, the branching points turn into first-order 
poles and this problem is absent (for details on these 
branch cuts see Ref. [i^). These branch cuts can be 


avoided if one takes the integration contour along the 
real time axis once the electron is in the continuum. 
More specifically, the time integration contour is taken 
first parallel to the imaginary time axis, up to the so- 
called “tunnel exit”, i.e., the point in space at which the 
electron tunnels out of the potential barrier, and then 
along the real time axis. This is the procedure taken 
by most groups when implementing Coulomb-corrected 
strong-held theories (see, e.g., Refs. [ssl. I^j. 

Furthermore, besides the SEA factor (I14|) and the tun¬ 
nel matrix element {pg(tg) + A{tg) \Hj{tg)\ tpo), the am¬ 
plitude now involves the stability of the trajec¬ 

tories. In the limit of vanishing binding potential, the 
usual SEA is recovered [H^. However, as we will see, 
this happens in a nontrivial fashion when the degeneracy 
breaking of trajectories is taken into account. 


III. RESULTS AND DISCUSSION 


In the results that follow, we use the monochromatic 
laser held 


E(t) = zEo sinwt. (24) 

For this type of held, A{t) oc cos(wt)z and given hnal 
momentum p, there are two solutions tg of Eq. (fT^ for 
the SEA per cycle of the laser held. In [i^, these solu¬ 
tions have been related to Orbits I and H, depending on 
whether the electron leaves in the same or in the opposite 
direction to the detector. In the results that follow, we 
will consider this classihcation and its extension to the 
Coulomb-corrected case. 

The initial states are taken as the ground state of 
Hydrogen, i.e., tpoir) = {r\tpo) = In this 

case, the tunnel matrix element in Eq. (1231) becomes re¬ 
lated to the atomic dipole moment and can be simpli- 
hed as (Ps(ts) + A(t.O | — r ■ E(C)| V'o) E{t)kz, where 
k = Ps{tg) -I- A{tg) [^. Unless otherwise stated, we will 
consider C = 1 in U(r). 


A. Coulomb-corrected saddle-point solutions and 
their physical implications 

In comparison with the SEA, the canonical momentum 
p of the electron is time-dependent according Eq. (EH) if 
the Coulomb interaction is incorporated. Therefore, in 
the CCSFA theory, the greatest challenge is to solve the 
saddle-point equations for the tunneling time t' and the 
canonical momentum pg for any given final momentum 
p. One should also bear in mind that, in experiments, 
the measured photoelectron spectra is a function of the 
final momentum. Therefore, if, for a given final momen¬ 
tum, the initial conditions for the corresponding electron 
trajectories could be obtained reversely, it would be eas¬ 
ier to understand how these trajectories were influenced 
by the Coulomb potential. 
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Polarization 



FIG. 1: (Color online) Illustration of the four types of CCSFA 
trajectories in the zx plane for electrons with fixed final mo¬ 
mentum p, computed using a linearly polarized square pulse 
of duration Tp — 15.25 cycles, intensity I — 2x lO^^W/cm^ 
and frequency lo = 0.057 a.u., and a Coulomb potential © 
with (7=1. The ionization potential has been taken as 
Ip = 0.5 a.u. The black dot denotes the position of the nu¬ 
cleus. The inset shows the region near the core. 


In order to simplify the calculation and isolate the 
main effects of the potential on the trajectories, we as¬ 
sume that the electron is ionized by tunneling from the 
time tg to = Rets at a fixed momentum Psi^s) and 
then moves to detector with the real time and coordi¬ 
nate according to the classical equations of motion (1^ 
and (E2I)- This is the most widely used ass ump tion for 
the contour, and has been employed in [35|, 1^, 1^ (for 
a review see 0)- Within this set of assumptions, it is 
a reasonable approximation to neglect the Coulomb po¬ 
tential in Eq. (fT51) . which thus reduces to Eq. (IT^ . The 
tunneling exit at the time Ir is given by 

zq = a{tf) — Re a{t ), (25) 

where a{t) = J* A(T)dr [23- 

Figure [T] depicts four types of the trajectories in the zx 
plane for an electron with a given final momentum p. For 
trajectories of type I, the tunneling exit zq > 0, and the 
electron moves directly towards the detector without re¬ 
turning to its parent core. For the type II and III trajec¬ 
tories, the tunneling exit z^ < 0, meaning that the initial 
motion carries the electron away from the detector before 
it turns around and ends up with the stipulated momen¬ 
tum p. A closer inspection shows that they are similar 
to Kepler hyperbolae to which a drift motion caused by 
the field is superimposed [H, Q • Trajectory types I and 
II are similar to the so-called “short” and “long” trajec¬ 
tories in the SFA theory. The type III is not found in the 
SFA and can be observed after the Coulomb potential is 
considered, which is consistent with earlier work [d^. 

As we show in Fig. [2Ka), the emergence of the new 
trajectory is directly related to the momentum non¬ 
conservation. It is instructive to start with the case of 





X (a.u.) X (a.u.) 




X (a.u.) X (a.u.) 


FIG. 2: (Color online) Emergence of the two distinct trajec¬ 
tories II and III. (a) For ionization parallel to the field, the 
trajectories are part of a torus, respecting the rotational in¬ 
variance of the system. In the SFA, this torus would contract 
to a point as momentum conversion then dictates that the 
initial momentum vanishes. At a finite tilt angle of the final 
momentum (here 2°), however, the torus splits into two dis¬ 
tinct solutions, II and III (red and green), (b) As the Coulomb 
interaction strength C is then reduced, trajectories II and III 
merge, and remain distinct from trajectory I (here the tilt an¬ 
gle of the final momentum is 1°). The remaining parameters 
are the same as in the previous figure. 


ionization parallel to the field, where the final momentum 
p is along the polarization direction. On first inspection, 
the type II and III trajectories are symmetric with re¬ 
spective to the polarization direction. As a matter of fact, 
trajectories II and III then degenerate into a torus, with a 
finite initial transverse momentum. In the SFA, the torus 
contracts onto a single trajectory with a vanishing initial 
transverse momentum, as dictated by momentum conser¬ 
vation. At a finite tilt angle, the torus splits in analogy 
to the Poincare-Birkhoff scenario in KAM theory, leaving 
two clearly distinct trajectories. As we change the effec- 
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X (a.u.) 


FIG. 3: (Color online) Behavior of type IV orbits for increas¬ 
ing photoelectron energy. The field and atomic parameters 
are the same as in the previous figures, and the energy is 
indicated on the upper right corner. A blow-up of the fig¬ 
ure shows the trajectories approaching the origin, which is 
marked by a black dot. The field and atomic parameters are 
given in Fig. [T] 


five Coulomb interaction strength C from (7=1 (Hy¬ 
drogen) to (7 = 0, the situation from the SFA is again 
recovered [panel (b) in Fig. [2]. 

It is noteworthy that our numerics uncover an addi¬ 
tional trajectory type, denoted as IV. Although the tun¬ 
nel exit points towards the detector, the electron is driven 
back to the core by the laser field, then goes around the 
core, and finally moves towards to the detector. With 
increasing photoelectron final momentum p, the short¬ 
est distance between the electron and the core decreases. 
This distance can be smaller than the tunnel exit. In this 
case, this type of trajectory corresponds to a rescattering 
event. This behavior is shown in Fig. [S] 

Since the contribution of events involving rescatter¬ 
ing to the low-energy photoelectron spectra is small, it 
can be safely neglected, so that we only need to consider 
the type I-III trajectories. However, it is encouraging 
to see that rescattering contributions already show up 
within a framework that is formally tailored for direct 
ATI, even though they eventually may require a more ac¬ 
curate treatment including rescattering form factors that 
account for the inherently diffractive, nonclassical nature 
of these events [H, [2^ . 


B. Photoelectron spectrum 

Based on the trajectories described in the previous sub¬ 
section we now study the photoelectron spectrum within 
the CCSFA theory and compare the result with the stan¬ 
dard SFA, while taking the ab initio TDSE calculation 
as a benchmark. The TDSE has been computed using 
the freely available software Qprop [tJ- The standard 
SFA is implemented according to Eq. (IT^ . Within the 
CCSFA, we calculate the stability of the trajectories nu¬ 
merically. In practice, instead of using 9ps(<)/9rs(ts) in 
Eq. (|2^ we employ 9ps(t)/9ps(is). The latter stability 



FIG. 4: (Golor online) Photoelectron spectra for ionization 
along the laser polarization direction, with the standard SFA, 
the Coulomb-corrected CCSFA, and from ab initio TDSE cal¬ 
culation. The blue solid curve is the envelope of the TDSE 
calculation. In order to perform a clearer comparison, the 
TDSE spectrum has been shifted upwards in around one or¬ 
der of magnitude. The laser-field and atomic parameters are 
given in Fig. [T] 


factor is of easier implementation, and can be obtained 
using a Legendre transformation in the transition am¬ 
plitude (1201) . Upon this transformation, the action will 
remain the same as long as the electron starts from the 
origin, which is one of the assumptions considered in this 
work. For details on Legendre transformations see, e.g., 


We use the tunnel exit approximation (l25)l to split the 
action into a part inside the barrier, 


STiP8^i^s,ts) = - H{p^,rs,T)dT 


(26) 


with the tunneling trajectory rs(t) = J* [p^(t)- f A(T)]dr 
and a part outside the barrier. 


ST\Ps,^s,ts) = [ dT[-p^(r) • rs(T) - iJ(p^,rs,r)], 

JtR 

(27) 

with the ionization trajectory determined as described in 
the previous subsection. Note that, in Eq. (1251) . the term 
in ps is vanishing as the momentum inside the barrier 
was taken to be constant. Furthermore, real variables 
outside the barrier will lead to real stability factors, so 
that phase differences in the continuum stem exclusively 
from the action. 

The resulting photoelectron spectra in the direction 
along the laser polarization are shown in Fig. |4l Since 
the SFA and CCSFA only accounts for the interference 
of the electrons ionized in one optical cycle, the spectra 
correspond to an envelope, without the sharp ATI peaks 
seed in the ab initio calculation. We therefore also show 
the envelope of the spectrum from the ab initio method. 
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Energy (a.u.) 



Energy (a.u.) Energy (a.u.) 


FIG. 5: (Color online) Photoelectron spectra from the 

CCSFA theory with different strengths C of the rescaled 
Coulomb potential given by Eq. ©• For C = 0 the spec¬ 
trum coincides with that of the SFA theory. The remaining 
laser-field and atomic parameters are given in Fig. [T] 


FIG. 6: (Color online) Amplitude of each trajectory type 
as a function of the photoelectron energy (a) within the SFA 
theory and (b) within the CCSFA theory. The laser-field and 
atomic parameters are given in Fig. [T] 


Clear interference structures are observed in all spectra. 
A closer inspection shows that the interference contrast 
in the spectra from the TDSE and the CCSFA theory is 
much weaker than that from the SFA theory. Moreover, 
the positions of the interference maxima in the CCSFA 
spectrum are in a better agreement with the TDSF result 
than the SFA. 

The mechanisms leading to these improvements are ex¬ 
plained in the next subsections. For reference, it is useful 
to inspect how the Coulomb corrections are established 
when one changes the effective interaction strength C, so 
that for C = 0 the SFA is recovered. Fig. [S] shows the 
photoelectron spectra for different values of C. 



Energy (a.u.) Energy (a.u.) 


1. Interference contrast 

First, we consider the effect of the Coulomb potential 
on the interference contrast. According to the discussion 
above, the type I-III trajectories are dominant for the 
electrons with low kinetic energy. For the photoelectron 
spectrum along the laser polarization, type II and III tra¬ 
jectories are symmetric with respect to the polarization 
direction and have the same phase and amplitude. There¬ 
fore, the interference pattern in the spectrum arises from 
the beating between type I trajectory and type II and III 
trajectories. In Fig. [HI we present the amplitude related 
to each orbit as a function of the photoelectron energy 
with and without Coulomb corrections, respectively. 

In the SFA theory, the amplitudes associated with the 
type I and the type II trajectories are the same. This 
holds because in the SFA the electron’s final momentum 
is solely determined by the ionization time, which for 
trajectories I and II are displaced by half a cycle. This 
means that the absolute values of the electric field, and 
hence the ionization probability, are the same. Therefore, 


FIG. 7; (Color online) Time of tunneling as a function of 
the photoelectron energy, (a) for the type I orbit and (b) 
for the type II orbit. The time is separated into two parts: 
real part and imaginary part. The black curves are for the 
SFA calculation and the red curves are from the CCSFA. The 
laser-field and atomic parameters are given in Fig. [T] 


in the SFA the interference contrast will be maximal. If 
the Coulomb potential is included, the amplitudes of the 
type I and II/III trajectories differ slightly. Furthermore, 
the joint amplitude of the type II and III trajectories 
exceeds that of type I significantly. All this leads to a 
much reduced contrast of the interference pattern [8l| . 

Phenomenological insight into the amplitude difference 
between the type I and type II/III trajectories can be ob¬ 
tained by considering the tunneling time and initial mo¬ 
mentum as a function of the photoelectron energy. The 
imaginary part Im [tg] can be interpreted as the time it 
takes the electron to tunnel through the potential barrier 
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FIG. 8: (Color online) Initial momentum for trajectories of 
type I (a) and II (b) as a function of the photoelectron energy. 
The momentum is separated into the parallel and perpendic¬ 
ular components. The blue curve in (b) denotes the total mo¬ 
mentum from the CCSFA theory. The laser-field and atomic 
parameters are given in Fig. [T] 


[g^. The larger Im is, the lower the ionization rate. 
Fig. [7] exhibits the time of tunneling as a function of the 
photoelectron energy for the type I orbit and the type II 
orbit, respectively. For the type I orbit, Im [tg] increases 
when the Coulomb potential is taken into account. In 
contrast, for type II and III trajectories, Im[ts] is smaller 
in the CCSFA than in the SFA. Therefore, if the Coulomb 
corrections are present, the amplitude of the type I will 
be smaller than that of type II. These features are consis¬ 
tent with the changes in Re[ts], which, for Orbit I, move 
towards the field crossing and, for Orbits II/III, is dis¬ 
placed towards the times for which the field amplitude is 
maximal. This implies that the effective potential barrier 
will be wider for the former orbits and narrower for the 
latter. 

According to the Ammosov-Delone-Krainov (ADK) 
theory 0 , these observations can be linked to the ini¬ 
tial momentum, with a large momentum translating into 
a reduced ionization rate. Fig.[5Ka) shows that for trajec¬ 
tories of type I, the initial momentum from CCSFA the¬ 
ory is indeed larger than that from the SFA calculation. 
This can be understood from the fact that the electron 
needs to compensate the deceleration in the Coulomb po¬ 
tential as it moves towards the detector. In contrast, for 
type II and III trajectories, the initial parallel momen¬ 
tum from the CCSFA theory is smaller than in the SFA. 
Although there is a nonvanishing perpendicular momen¬ 
tum from the CCSFA theory, the total initial momentum 
[see blue curve in Fig. [Hb)] is still lower than that from 
the SFA theory. This indicates that the electron accel¬ 
erates significantly due to the interplay of the Coulomb 
potential and the laser field as it passes near the core. 
We conclude that the amplitude difference between the 
different types of trajectories is generally consistent with 



Energy (a.u.) 


FIG. 9: (Color online) Phase difference A<I> = <I>i — "F// (in 
units of tt) between trajectories of type I and II as a function 
of the photoelectron energy. The arrows denote the positions 
of the interference maxima, where AF = 2k'n, k = 1 , 2 ,.... 
The laser-field and atomic parameters are given in Fig. [T] 

the phenomenological picture of the ADK theory. 

Physically, the above-mentioned behavior can be at¬ 
tributed to the fact that the Coulomb potential deceler¬ 
ates the electron for Orbit I, which hinders ionization. In 
contrast, for orbits II and III, the Coulomb potential ac¬ 
celerates the electron. Hence, the electron acquires an ad¬ 
ditional pull, and may escape moving along laser-dressed 
Kepler hyperbolae. 


2. Positions of interference maxima 

We turn to the positions of interference maxima in the 
spectra. Fig. 2] shows that the positions of the interfer¬ 
ence maxima in the spectra from the CCSFA theory are 
shifted when compared with the SFA. Since the positions 
of interference maxima are determined by the phase dif¬ 
ference between different types of trajectories we study 
how this is affected by the Coulomb potential. In Fig. [HI 
we show the phase difference between the type I and type 
II as a function of the photoelectron energy. After con¬ 
sidering the Coulomb potential, the dynamical phase dif¬ 
ference from the CCSFA becomes smaller than that from 
the SFA. This can be traced back to the fact that the 
type II trajectory accumulates a larger negative phase 
contribution from the Coulomb potential as it passes by 
the core. Overall, this reduces the phase difference to 
trajectory I. A similar analysis has been employed in our 
previous publication [2^, in the context of ATI with el- 
liptically polarized fields. 

Note that the reduction of the phase difference ap¬ 
proaches 27r, so that neglecting multiples of 27r one can 
also interpret the large shift towards lower energies as 
a small shift towards larger energies. This ambiguity is 
resolved when we consider the effect of the continuously 
rescaled Coulomb potential in Fig. [51 which shows how 
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FIG. 10: (Color online) Photoelectron spectra from the 
CCSFA theory and the SFA theory, respectively. The red 
curve is the spectrum from the CCSFA, but without the sub- 
barrier Coulomb term I = — V[r{r)]dr. The laser-field 
and atomic parameters are given in Fig. [T] 


a large shift towards smaller energies is established as C 
increases. These results are also consistent with the re¬ 
cent TDSE simulations in [^, and with the outcome of 
similar Coulomb corrected approaches [H, . 


3. Overall ionization amplitude 


Finally, we turn to the overall ionization amplitude 
of the spectra. As we can see from Fig. 01 the ion¬ 
ization amplitude from the Coulomb-corrected theory is 
much larger than that from the SFA, and much more 
in line with the TDSF result. Enhanced tunnel ioniza¬ 
tion is a well-known effect caused by Coulomb corrections 
to the effective potential barrier. It was first predicted 
in and subsequently observed in several Coulomb- 
corrected strong-held calculations [ill, [t^ . Furthermore, 
in the 1980s the Coulomb-induced orders-in-magnitude 
enhancement of tunnel ionization rates of atoms and pos¬ 
itive ions was well documented in experiments [o^ . In¬ 
tuitively, it can be understood that compared to the SFA 
theory, the barrier for a Coulomb potential is smoother 
and lower, making it easier for the electron to tunnel 
through. Within the CCSFA, this enhancement of the 
ionization amplitude mainly originates from Coulomb 

term / = — V'[r(T)]dr in the sub-barrier action (E51) . 
Fig. (TU] shows the result if this term is neglected. The 
overall magnitude of the spectrum is then comparable 
to the SFA. In the CCSFA this term contributes with a 
negative imaginary part ImJ < 0, and thus increases the 
ionization amplitude. 


IV. CONCLUSIONS 


In this work, we develop and use a path-integral for¬ 
mulation to assess the influence of the residual Coulomb 
potential in above-threshold ionization. We focus on 
the direct transition amplitude, in which hard collisions 
with the core are not incorporated. Overall, the pho¬ 
toelectron spectra obtained with the Coulomb-corrected 
method presented in this paper exhibits an excellent 
agreement with the ab initio solution of the time de¬ 
pendent Schrbdinger equation, which is far superior to 
that encountered for the plain strong-held approximation 
(SFA). This is especially true for the interference sub¬ 
structure in the spectra. We also perform a systematic 
analysis of how the Coulomb potential modihes the or¬ 
bits along which the electron may leave its parent ion and 
reach the detector. We compare our results with those of 
the SFA, and make an assessment of how, in the limit 
of vanishing Coulomb coupling, the SFA is recovered. 
The present formulation is closely related to the concept 
of quantum orbits widely employed in semi-analytical 
strong-held approaches with and without Coulomb cor¬ 
rections. 

We have built upon the existing knowledge that the 
Coulomb potential introduces a richer topology in the 
electron motion [d^, with four distinct sets of orbits, and 
have related these orbits to those in the SFA in a more 
systematic way. Throughout, we have employed the same 
classihcation as in [dj, which specify these four sets as 
Orbits type I to IV. In particular, we have found that, for 
electron emission along the polarization axis, due to the 
rotational symmetry with regard to the field axis. Orbits 
II and III will be located on a torus. This torus will con¬ 
tract for decreasing Coulomb coupling, until it becomes 
a point. Physically, this means that, for the SFA, Orbits 
type II and III will merge into a single, degenerate orbit 
if the final electron momentum is parallel to the laser- 
field polarization. For non-vanishing emission angle, the 
above-mentioned torus will break down, and there will 
be two discrete solutions. This behavior indicates that, 
strictly speaking. Orbits type II and III should not be 
treated independently in an asymptotic expansion when 
computing photoelectron spectra and momentum distri¬ 
butions. Indeed, a rigorous treatment would require solv¬ 
ing the integral around the manifold exactly for a final 
momentum along the axis, and a uniform approximation 
for non-vanishing emission angle \l‘A. l73l . This strongly 
suggests that the cusps observed in [^,13 close to the 
so-called ATI low-energy structure are related to this ef¬ 
fect. It is indeed noteworthy that a very good agreement 
between the TDSE and the Coulomb-corrected SFA in 
[sM was obtained throughout, except in this region 
(see also the discussion of this cusp in the review [671]). 
We have however not studied the above-mentioned cusp 
systematically. 

Furthermore, our results indicate that, if the Coulomb 
potential is accounted for, the concepts of “direct” and 
“rescattered” electrons are not very clear-cut. These con¬ 
cepts are very clear in the SFA, as there are either hard 
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collisions with the core, or no collisions at all. If the 
Coulomb corrections are present, however, the Coulomb 
potential strongly deflects Orbits type IV. These orbits 
go around the core, and there is a marked decrease in the 
electron’s shortest distance from the origin as the pho¬ 
toelectron momentum increases. For high enough mo¬ 
mentum, this distance is located in a region in which the 
binding potential is dominant. This behavior could be 
interpreted as a type of recollision, which is absent in the 
SFA. The amplitude associated with this type of trajec¬ 
tories is however very small and hence not relevant to 
the computation of ATI spectra in the parameter range 
of interest. 

In addition to that, we have investigated the influence 
of the Coulomb potential on the ATI spectra, with em¬ 
phasis on the interference contrast and position of the 
maxima. This influence has been traced back to partic¬ 
ular sets of orbits. First, the contrast in the interfer¬ 
ence structure decreases, in comparison with the SFA. 
This happens because, in the SFA, Orbits I and II are 
equivalent, and displaced by half a cycle, while, if the 
Coulomb potential is included, this no longer holds. In 
fact, the Coulomb potential will decelerate the electron 
if it reaches the continuum along Orbit I, and will accel¬ 
erate the electron if it is ionized along Orbits II and III. 
This will lead to an increase in the amplitudes associated 
with Orbits II and III, and to a decrease in the ampli¬ 
tude related to Orbit I. Furthermore, there is the joint 
effect of Orbits II and III, which will weaken the fringes. 
Recently, the influence of Orbit III on interference effects 
has also been investigated in a different context, namely 
side-lobes in ATI electron momentum distributions, and 
it has been found to be significant [t^ . 

The suppression of Orbit I and the enhancement of Or¬ 
bits II and III has been confirmed by a systematic anal¬ 
ysis of the initial momenta and ionization times. For Or¬ 
bit I, the initial momentum increases when the Coulomb 
potential is considered, while, for type II/III orbits, the 
initial momentum decreases. Physically, this means that 
the Coulomb potential hinders ionization along Orbit I, 
as the electron will require a larger momentum to escape. 
For Orbits II/III, the Coulomb potential accelerates the 
electron after the tunneling ionization, so that a lower 
escape momentum is required. An increase in the ini¬ 
tial momentum for the type I orbits also implies that the 
electron ionization time has moved away from the field 
maximum towards the field crossing. This means that 
the effective potential barrier through which it must tun¬ 
nel will widen. Hence, there is also an increase in Im[ts]. 
In contrast, for Orbits type II/III, the tunneling time 
moves to the crest of the laser field and thus the effective 


potential barrier becomes narrower. These observations 
are consistent with the changes in the real parts Re[ts] of 
these times, as shown in Fig.jT] 

Similarly to the results reported in [H, , we also 

find that there is a phase shift towards lower energies in 
the interference maxima. In our model, this phase dif¬ 
ference occurs due to Coulomb effects in the continuum 
propagation, while sub-barrier corrections mainly influ¬ 
ence the overall yield. In contrast, in [11,14^, this phase 
difference is attributed to sub-barrier corrections instead. 
While the contour taken by us and the assumption that 
all variables are real outside the barrier are also employed 
in [H, the term Ps(r) • rs(r) is absent in their ac¬ 
tion. Eq. (EU shows that this term is proportional to 
the gradient of the binding potential. Its value is small 
for Orbit I, which moves towards the detector directly, 
while it is large for Orbits II and III, which are deflected 
by the core before reaching the detector. We have indeed 
verified that the phase from this term plays an important 
role in our formulation. Indeed, if this term is removed 
from the action, there is significant deviation between 
the TDSE and CCSFA results. The stability factors em¬ 
ployed here are also different from those in jss], , 

but they influence mainly the contrast and not the posi¬ 
tion of the maxima. They are, however, very important 
for a quantitative agreement with the full TDSE com¬ 
putations. This term is also absent in [b^, in which 
the Eikonal-Volkov approximation is employed and the 
phase differences are obtained along propagation by us¬ 
ing a complex intermediate coordinate. Since, however, 
in circularly polarized light is used, it is expected 
to be vanishingly small as the electron never returns to 
the core after tunneling ionization. We expect that the 
present analysis will contribute to a better understanding 
of cusps and the ATI low-energy structure in the future. 
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